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Abstract 

Expressing Weierstrass type infinite products in terms of Stieltjes inte- 
grals is discussed. The asymptotic behavior of particular types of infinite 
products is compared against the asymptotic behavior of the entire func- 
tion £(s), well-known in Riemann zeta function theory. An approximate 
formula for the distribution of the non-trivial roots of Riemann's zeta 
function is obtained. 
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1 Introduction. 

In the theory of Riemann's zeta function the entire function 

£( s )4r^V* s ( s -i) C ( s ) (i.i) 



2 V2, 

plays a fundamental role. As demonstrated in this theory, the zeros of £(s) 
occur only in the critical strip < 5R(s) < 1, and given the so called Riemann 
hypothesis is true, these zeros will all be located on the vertical axis 3?(s) = \. 
Theorems developed by Hadamard in 1893 allow to demonstrate that £(s) may 
be represented by an infinite Weierstrass product (see e.g. Lang, p. 386): 

e(s)=e(o)n( i - £ ) (i-2) 

p \ p / 

in which the set p represents the zeros of £(s), which are all located within 
the critical strip (see e.g. Edwards, p. 18). An important point to note in the 
product l|1.2p is that no exponential convergence factor is required per individual 
term. £(s) is analytic in the entire s plane and is an even function of s — |, i.e. 
£(s — |) = £(| — s), and its roots occur in pairs, i.e. if the complex number p 
is a root, so is 1 — p. 

In (|1.2|) the terms in the product are to be taken in pairs, i.e. p together with 
I- p. 
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2 Types of infinite products. 

In what follows we will consider two types of infinite products which are a 
generalisation of (|1.2p : 

1. Pure axial products in which all zeros lie on the axis 5ft(s) = i. 

2. Strip type products in which some zeros also lie outside the same axis as 
in 1 but still within the critical strip. 

The axial products are a subset of the strip product set. £(s) belongs at least 
to the strip product set and if the Riemann hypothesis were true, £ (s) would be 
an axial product. For both 1 and 2 we exclude roots at s = \. 
Below, instead of the complex variable s = a + it usual in analytic number 
theory, we will use: 

1 

z = s— — = x + it 

This brings more symmetry into the considerations. For both axial and strip 
types we will consider products F(z) analytic in the entire z plane (entire func- 
tions) which are in addition even, i.e. 

F{z) = F(-z) 
which implies F[x) — F(—x) and F(it) = F(—it). 

We also assume F{x) as well as F(it) to be real and the product expansion for 
F(z) requiring no exponential convergence factors. We will therefore consider: 

1. Axial products: their zeros are purely imaginary and occurring in pairs 
±z;with: 

zi = ih h>0 (2.1) 

The product 

oo 2 

F( Z ) = F(0)n(l+p) (2-2) 
1=1 1 

with F(0) real, represents an even entire function in the z plane. As well 
known 112.21) converges if 

(2-3) 



1=1 



converges 



2. Strip products: their off-axis zeros need to occur in quads, i.e. z m , — z m , 
~z m ~,-~z m ~, all four values lying within the critical strip. We will use the 
representations: 

- — j/3m 



iq m e-^ m (2.4) 



2 



■ i, 



in which we assume q m to be real and > \ 
The associated product then equals: 

°«-^n('~)K)(i~)H) 

with G(0) real. 

As for axial products we assume X)m=i to converge. Strip products then 
consist of the combined product of expressions (|2.2p and j|2.5|) • In what follows 
we will first consider axial products. Later we will look at strip products and 
multiple roots. 

3 Axial products. 

By taking the logarithm of H2.2j) we split the product into a sum: 

oo 2 

f (z) =lnF (z) =lnF (0) + ^ln(l+^) (3.1) 

1=1 1 

In F(z) is here defined by: 

f z =JnF + / -F^dw 
Jo F ( w ) 

in the open semiplane x > 0, via a path of integration not crossing the axis 
x = 0. This way f{z) is defined unambiguously and the singularities occurring 
at z — ±iki are avoided. The fact that F(0) is taken at the point z = 0, i.e. 
on the axis x = is no problem for the definition of f(z). In conjunction with 
the distribution of the zeros z = ±iki we now introduce a non-decreasing step 
function for k > 0: 

¥>(k)= E 1 ( 3 - 2 ) 

ki<k 

In words: for fc; < k < fej + 1, y(fc) equals the number of roots less than or equal 
to ki. The way we have defined our infinite products, ip(k) always equals zero 
for < k < k% and then jumps to the value 1. In order for (|3.ip to converge, 
the infinite sum T? should converge (see l|2.3p ). 

This implies also that there is no accumulation point for the ki for any finite 
value of k. We also ignore multiple roots for the moment. (|3.ip may now be 
expressed as a Stieltjes integral: 
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f(z)=ZnF(0) + / ln(l+—)d<p(k) 



(3.3) 



Partial integration yields: 



f (z) =lnF (0) +y (k) In ( 1+^ 



Jo 



equals 



For the second term on the right, in the lower limit (p (k) In 
zero since <p(k) = for < k < k%. 

Consider now the upper limit (k — > oo) of the same term. For a given value of 
z and k^ oo: 



z 2\ z 2 



As long as for k — * oo, y>(fc) < C k a with C a positive constant and a < 2, the 
term (^(k^n ( 1 + p-J goes to zero as well. This is a constraint we put on tp(k) 



and is in fact the equivalent of (|2.3p Under these assumptions: 

r,2 i 



f (z) =lnF (0) 



p(k)d 



f (z) =ZnF (0) + 2z 2 



ln(l+ 
V?(k)dk 



k 2 



/o k(k 2 +z 2 ) 

4 Axial product example: Cosh z 



(3.4) 



A straightforward example is the function F(z) — Cosh z, which for z = it 
becomes F(it) — cos t. 

Cosh z is an entire function and its zeros occur in pairs along the imaginary 
axis at equidistant points: 



ki= ±(2i-i)- a = i,2,...) 



and as well known: 



Also 



Coshz=IJ(l 



4z 2 



\X t 2 (21-1) 



(4.1) 



Coshz 



e z +e" 



(4.2) 



We will now show how to recover l|4.2p from (|4.ip . through the use of p.4p . In 
since F(0) in (f47T|) equals 1, In F(0)=0, so in this case becomes: 

2 r v( k ) dk 



f (z) = 2z- 



o k(k 2 +z 2 ) 



with 



k>(21-l)f 



(4.3) 
(4.4) 
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We split l|4.4p up into a continuous part: 



and a sawtooth type odd periodic function tfi2{k) with period it which over the 
interval — 5 <fc < ^ equals: 

ip 2 (k) = -- 

7T 

For k > -| ^2 (fc) just repeats itself. 

For our purposes we only need to consider <fi2{k) in the range < k < oo . 

¥>(k)=-+pa(k) (4.5) 



(f2(k) may be expanded in a Fourier series: 

oo 

^ 2 ( k )=]T(-ir 



sin 2nk 



(4.6) 



From 03 



f (z) = 2z 2 / 
Jo 



dk 



k 

-+^ 2 (k) 

7T 



k(k 2 +Z 2 ) 

The Dirichlet conditions apply to (p2 (k) , so we are allowed to integrate term by 
term for ip2(k). 

dk ^2, (-1)" f°° sin(2nk) dk 



rn 2z 2 



o (k 2 +z 2 ) 



k(k 2 +z 2 ) 



(4.7) 



The integral J °° y(y£ + a%) may be obtained via residue techniques and equals 



Through appropriate change in variables: 
2z 2 ^2, (_l) n f°° si n2nk dk ^ (-l) n 



; z 2 - (_!)° ro 



k(k 2 +z 2 ^ 



n=l 



(1- 



oc / ., \ n oo / . \n+l 

.) = y L1L+ y izli — e - 2nz 



n=l 



Both series are conditionally convergent and equal respectively — ln2 
and ln{\ + er 2z ). (|4.7p then becomes: 
f{z) = z-ln2 + ln(l + e~ 2z ). 

Through exponentiation and analytic continuation to the entire z plane: 

pZip — z 



(4- 
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5 The transform from <p(k) to f(z). 



For a given ip(k), the zeros of f(z) are unaffected by the assumed value of F(0). 
Conveniently in this section we simply put F(0) = 1. (|3,4p then reduces to: 

'<■>-* fi^j ^ 

To the step function <p(k) as defined by l|3.2p corresponds a specific f{z) as given 
by l|5.ip and a specific F(z) = e^ z ' as well. 

Assuming F(z) constitutes an axial product as defined in this paper, ip(k) may 
be recovered unambiguously from F(z) by taking half of the integral 

1 I F '^-dz 



2vri J Ck F(z) 

where Ck is a contour dependent on k in the sense that it encircles the imaginary 
axis counterclockwise, crossing this axis ones at each of the points z = ik and 
z = —ik, e.g. a circle with radius k and centered at z = 0. 

We should in this context keep in mind that although f(z) = InF(z) obtained 
through (|5.ip is restricted to the semiplane x > 0, exponentiation to F(z) and 
analytic continuation yields an analytic function F(z) over the entire plane. 



6 Transform table. 

In table 1 are listed the results of the transform (|5.ip for a number of different 
functions <p(k). With the exception of a possible initial discontinuity, all these 
functions are continuous for the rest of the interval up to oo. Although tp(k) 
as defined in l13.2[) is non-decreasing, some of the functions in table 1 are non- 
increasing or decreasing. By themselves such tp(k) do not make sense. Added to 
a non-decreasing (p(k) however, they do make sense. Functions similar to f2(k) 
in I|4.5J1 , which have discontinuities over the entire range of k are not included in 
the table. Although such contributions to the total (p(k) are essential in order 
to create zeros for F(z), for the asymptotic behaviour of F{z) however they just 
add a multiplicative constant. This is well illustrated by the example of Coshz 
in section 4: For — > oo, Cosh z ~ whereas for the zeros of Cosh z, the 
contribution of (f2(k) is essential. Without it, %- would not have any zeros. In 
addition it is essential for making F{z) an even function. The ip(k) listed in 
table 1 are more than sufficient for our further needs. The function u(k — a) in 
table 1 represents the unit step function, i.e. u(k — a) = for k < a, while it 
equals 1 for k > a. 
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Table 1 |[ 








<p(k) 




/(*) 


1 


u(k — a) 


a>0 


Zn (l + fr) = 2/nz - 2/na + (z~ 2 ) 


2 


ku(k — a) 


a>0 


7tz - 2a+ (z~ 2 ) 


3 


lnku{k — a) 


| z |> 1; a>0 


(lUZ) + 12 2 2^ m= l m S z 2m 


4 


klnku(k — a) 


| z |> 1; a>0 


7rzZnz — 2a(/na — 1) + (z~ 2 ) 


5 


/nk.u(k — a) 


| z |> 1; a>0 


2 


, (-l) m+1 ?tlnz 


k 


Z J m=l (2m-l) 2 z J » z 


6 


kVk 




TTV^Z-y/z 


7 


k\fkln k 




7r ^ z (2Znz+7r) 


8 


u{k — a) 
k 


a>0 


2_i + 2 r a rctg(^)l 

a. z z o V z y 


9 


u{k — a) 
k 2 


a>0 


a 2 7. 2 \ a 2 ) 



7 Strip products 

In strip products some zeros lie off-axis. The total product equals H{z) — 
F(z)G(z) in which F(z) contains the axial zeros and G{z) the off-axis zeros. 
Let g(z) = In G(z). 
From (2HJ): 

g(zMnG(0) + f>+ ^5^ + ^) 

m=l q ™ qm 

For strip products we need to modify the definition of InG(z) slightly since zeros 
of G(z) occur for — \ < x < \. We therefore consider the definition of G(z) 
only for the halfplane x > \ . The formula used in section 2 for the logarithmic 

function g (z) —lnG(0)+ J °° %^ydiv is still valid if the path of integration starts 

at z = 0, follows the real axis up to z = i, and from there onwards to a point z 
for which x > \ 

For the asymptotic behavior of g(z) this restriction poses no problem. Let's 
assume a quad of zeros with parameters q m and /? m . Let C(q m ) be the circle 
with radius q m and centered at z — 0. Let P be the intersection of C(q m ) with 
the vertical z = ^ + it, with i > 0. 

Let /?m be the angle between the vertical z = it and the radius from z = to 
P. Then: 

(7<1) 

The first zero in the quad (see (2.4)), i.e. !g m e" 1 '' m will be located on C(q m ) 
and the associated angle /3 m will be less than 0m . Therefore: 

cos m >cos 0m (7.2) 
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jZTIJ) implies 



(f772|) and (f773|) imply: 



COS /?Af > 1 



1 > cos(3 m >1 



(7.3) 



or: 



cos (3 m =1 — — ^- wii/i < # m < 1 



(7.4) 



7]4| implies: 



In 1 



2z 2 cos/3 n 



-r = 2ln[ 1+ — 



-In 



4c4(l+f ) 



(7.5) 



For 5R(z) — » oo, we assume 

3f?(z 2 ) > Qi.e. |arg(z)| < 



(7.6) 



7.6p implies In 



1 - 



will go to zero for $t(z) — > oo. 



Therefore all finite sums of such terms will equally go to zero as 0(z 
Let's now take the case of an infinite sum of q' m s. 



Under constraint 117.61): 



It is not difficult to prove that in this case and for q m > \ 



< 



y c 



In 



1 



4^(1+^") 



<K with K some positive constant. 



For a strip product In H{z) = In H(0) + f(z) + g(z) 
with: 



f(z)= X>(l+p) 



1=1 



g(z) = ^2/n(l+^)+ J> 

1 4m „_i 



z 2 n 



4q m (l+^) 



As discussed above, the contribution of the second term of g(z) will be 0(z 2 ) 
for a finite number of terms and at most O(l) for an infinite number of terms. 



S 



The first term of g(z) is simply the equivalent of a double axial product term, 
i.e. each q m may be treated as a jump of size 2 in the <p(k) function. In fact, 
the same argument applies to multiple roots. Multiple roots will simply cause 
larger jumps in ip(k). 

As an example of multiple roots we consider an exercise similar to the one for 
Coshz (section 4). Here we consider <p (k) = ^+(fi2 (k) where for — < k < 
f2 (k) = — - and for ^ < k, (f2 (k) repeats itself with period Nn. N being a 
positive integer. The resulting ip(k) has on the average over a given interval 
< k < K, the same number of zeros as for Cosh z, but now we bunch N 
zeros together, i.e. the size of the jumps equals TV instead of 1. Each zero has 
multiplicity N. 

A similar calculation then yields: 

f (z) = z - N/n2 + NM(l+c- 2 t) 



F(z) = 



,f(z). 



Cosh 



For z — > oo: 



[Cosh(£) 



1 



Cosh z 2 N -! 

So for this example, bunching the zeros together leads to an asymptotic reduc- 
tion of F(z) equal to . 



8 Asymptotic behavior of In £(z). 

Note: It would be misleading to write £(z+^) for £(s) after substitution s = 
z+\. Indeed £(z) is an even function of z, not of z + ^. 

e(s)=^r(|)7r-is(s-l)C(s) 
For In r(|) we use Stirling's formula (see Lang, p. 423): 

InT (A) = ( A -0 lnA - A+^n (2tt) -w (A) (8.1) 

For A we will need to plug in A =|=f + |; 
The function w(A) equals: 

,^ /' 00 P(A)dA , N 

w<A »=/ -sir m> 

in which P (A)is a sawtooth type odd periodic function with period 1. In the 

interval < A < 1 P(X) equals — \ + A. After that it repeats itself. 

w(A)is analytic in the complex A plane from which the origin and the negative 

real axis is removed. 

w(A)m&y be written in various forms: 
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oo oo 



1 



m + 1 2m 



w(A)=-EE 

r=l m=2 

or more condensed (see Freitag-Busam, p204): 

oo 

w(A) = -£ 



1 



(A + r) 1 



k=0 



A + k 4) ln ( 1+ ATk 



or using the generalized or Hurwitz zeta function, 

1 oo _ 

C g (2m,A+-)=^(A+-+k) 

k=0 

C g (2m, A+i) 



w(A) = V , 
v ' ^ (2m + 

m— 1 v 

It is easy to demonstrate that for !ft(z) > 8^: 



1)2 211 



w 



z 1 

2 + 4 



<l x=» W 



(gljlyields: 

<" r (f+i 

Further: 



z /z 
2 V2 



4 V2 



-i 



4 te 1+ 2^ 



(8.3) 



z 1 1 , , /z l x 
.-- i+ -^( 27 r)-w -+- 



(8.4) 



ln£ (z) = -Zn2+Znr( -+-)- ^-+-J Zmr+ln ( ^--J +Zn£ (^z+-J (8.5) 

Due to the presence of Zn£ (z + 5) and In (z 2 — j) we define m£(z) only for the 
semiplane x > \ , similar to what is done in section 7, but excluding z = \ itself. 
Considering that for x > 10 we have \ln((z + l/2)\<-^ and putting some 
bounds on series expansions of certain logarithmic terms, we obtain: 

, , Z f Z *\ Z 7, 1, /7T\ 2z?i 

/n£(z)=-/n — h- Znz+-^n - H (8.6) 

^ W 2 V27r/ 2 4 4 12/ z V y 

with < l,$R(z) > 10 

At this point one should not get the impression that ((s) is not significant since 
its contribution to the asymptotic form of In £(z) disappears in a term 
Indeed, ((s) is an essential component in I|1.1J1 to create all the zeros of £(s) and 
making it an even function. 



9 Distribution of the roots for £(z). 

The intention now is to find an approximate formula for the number of roots up 
to a certain k value, based on l|8.6|l . To keep things simple, in what follows we 
will consider real values of z. This will not constrain our conclusions. 
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I|8.6p contains three terms which keep increasing without bound in absolute 
value with increasing z: 

Ti = —In ( — ] ; T-2 — ; T-i = —In z 

2 \2tt) ' 2' J 4 

In addition there is a constant and a vanishing term for z — > oo: 

If is the actual step function representing the distribution of the roots of 
£(z), we can split (p(k) up into a continuous curve (f>(k) and a discontinuous 
sawtooth type function Cl(k): 

p(k)=^(k)+n(k) (9.1) 

d = ~ d dic^ f° r a l mos t all k except for k — ki values corresponding to roots 
of £(z). At those ki, jumps will occur in f2(fc) equal to the multiplicity of the 
root. We will first concentrate on the terms T\,T2,T^. 

The idea is to identify for each of these terms a corresponding term in table 
1. We are of course allowed to multiply the <p(k) in table 1 with appropriate 
coefficients in order to obtain the required correspondence with (|8.6j) . 
To: 



k z a 
—— Ink.u (k — a) corresponds 77^ nz (Ina, — 1) +0(z~ 



k z a 

In (2tt) .u (k — a) corresponds — ln2n-\ — Zn27r+0(z ) 

2ir 2 tt 



k z a 

— — u(k— a) corresponds — -H hO(z ) 

2tt 2 tt 

7 7 7 

— u(k — a) corresponds — Inz — Zna+Ofz -2 ) 

8 4 4 

In total to: 



_2tt 




k 


7" 


v 2vr) 


~2tT 


h 8. 



.u (k — a) corresponds 

z / z \ z 7 a 7 
t:^ ( tt - I — 7: + t'"- z H — [2— ^na+Zn27rl ~-ln& + Ofz 
2 V27T/ 2 4 tt L J 4 v 

For k > a we may drop u(A; — a). To shorten things we call: 



-2\ 
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T <«4 n (£)-| + 4 inB (9 - 3) 

a 7 

T 5 (z)=T 4 (z)+- [2-lneL+ln2ir]--ln&+ 0(z~ 2 ) (9.4) 

7T 4 

We select the value of a such that we start with </>(a)=0. This happens for 
k = a = 9.6769... 
T§ (z) then becomes: 

T 5 (z) = T 4 (z) + 0, 8582 • • • + O (z~ 2 ) (9.5) 
We now compare T5 (z) with ln£(z), keeping in mind that 



00 / 9 
z z 



1=1 



(z) =Jn£ (0) +ln JJ ( 1+p ) (9.6) 



In £(0) may be computed from ITT]): Zn£(0) = -0,69892... 

From (EH) and J9j|: 

00 / 2 \ 1 

/n n ( i+ ^ ) - t s ( z ) = i /n (|) (°) -°> 8582 + ° ( z ~^ 



1=1 



In Y[ 1 J -T s (z) = 0, 0464 + O (z^ 1 ) (9.7) 
1=1 ^ 1 ' 



The remaining constant 0,0464 in l|9.7p must be due to fl(k), i.e. 

n(k)dk 



lim^ 00 2 /°° -£f^£r should equal 0,0464. 

fe( 1H 13- ) 



10 T/ie discontinuous function Q(k). 

In line with earlier definitions (see(3.2) . l|9.1|) ) we may put Cl(k) = for ^ 

k ^ a in which a is the point where <p(k) starts off from zero. According to the 

foregoing, a < ki where ki corresponds to the first zero of £(z). 

Let in general the transform 1 from the k domain to the z domain be understood 

as: 

ft (k) dk 



Let: 



Let also: 



2>(k)] = 2z 2 I j^-^y (10.1) 



limiwoo^ J ft(k)dk=ft (10.2) 
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Q (k) =fii (k) +fi (k) with n (k) =fi u (k - a) 



(10.3) 



T [fi (k)] =1 [fii (k)] +T [fi u (k - a)] = % [fii (k)] +O ln ^1+^ J (10.4) 

The existence of the logarithmic term in (|10.4|) as part of ln£(z) is in contradic- 
tion with (|9.7p where no logarithmic term is left. So O must equal zero. 
Therefore by definition: 

lim T ^oo / Oi(k)dk = (10.5) 

J a 

If fii (k)were a smooth function, it could after one or more zero-line crossings 
remain permanently above or below this line and approach it asymptotically 
and still meet (|10.5|1 . 

But Qi (k)is a discontinuous function with jumps of at least height one. There- 
fore Qi (k) will cross the zero-line an infinite number of times. If we now look 
back at the situation in terms of (p(k) and 4>(k) (see l|9.ip ). ip{k) will cross <j>{k) 
alternatively vertically and horizontally. 

At the horizontal crossings the number of zeros of £(z) from k — to the point 
of crossing k a is equal to (f){k a ). If two consecutive horizontal crossings occur at 
k a and kj,, the number of zeros in between will equal: 

k a +k b 



0(k b )-0(k a )~^^ln 



47T 



with A (a, b) =k b — k a . 

At the horizontal crossings the total number of zeros from up to T will equal: 

T ( T \ T 7 

N ^ = ^{^)-^ + s (10 " 6) 

In (|10.6p the symbol T has been used instead of k, as customary in analytic 
number theory. 

The reader is invited to go through an exercise in graphics. Draw the following 
curves in a diagram: 

1. <f)(k) starting with k = 9.67... and up to e.g. k = 55. 

2. </?i(fc): The tp(k) corresponding to £(z) based on the actual zeros available 
in the literature (e.g. Edwards, p. 96) up to e.g. the tenth root. <pi(k) 
is a step function with vertical jumps at the zeros of £(z). The horizontal 
portions occur for vertical values equal to n=0,l,2,... 

3. <fi2(k): A step function similar to tpi(k) with horizontal portions at the 
same level n=0,l,2,... as <pi(k), while the vertical jumps occur at those 
horizontal values k for which <j>(k) reaches a value n + |. 

It is remarkable how well tf2{k) approximates (fi(k). 
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